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In this paper solutions are presented for the determination 
of transient thermal stress distributions, in holiew circular 
eylinders of elastic material, in which the temperature is a 
function of the radial coordinate as well as time. Indicated 
Solutions are given for an externally insulated cylinder where the 
surface temporature is a known function of time, and for the case 
where the temperature-time relations are known for a heating fluid 
within the cylinder. Wuamerical results have been included for the 
specific case of a hollow steel cylinder subjected to the flow of 
saturated steam. For the particular set of conditions considered 
the tancential and axial stresses were found to be @ maximum at 
the heated surface at approximately five seconds after the start — 
of heating and attained a magnitude of 52,500 pounds per square 
inch. Stresses of this magnitude if often repeated could cause 
fatigue failure and result in rupture of a line carrying any 
appreciable pressure. It is conceivable that for more extreme 
conditions the surface stresses might reach magnitudes high enough 
to cause & permanent set in the material and residual stresses 
would be present when the cylinder was again cooled to the initial 
conditions. Stresses of these magnitudes warrant consideration 
in the design of steam lines that may be subjected to repeated 


heating cycles. 
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INTRODUCTION 


The thermal stress distribution in hollow circular cylinders 
cansed by & temperature gradient in the radial direction only is 
well known. The case of steady state heat flow in such cylinders 
is treated in several textbooks. The problem of etrese distri~ 
bution for the case in which the temperature is a function of both 
the radial coordinate and time lends itself to solution by the same 
stress equations, providing that the temperature-time relations are 
known. Analytic solutions for transient temperature distributions 
are given by many authors for the case where, from an initial uni~ 
form temperature, one or both of the surfaces are subjected to 
linear temperature gradients, or are suddenly raised to a civen 
value and thus maintained. The process for solution of such prob< 
lems is to substitute inte the already evailable stress equations, 
expressions for the temperature aga a function of time and radius 
for the civen condition of heat flew. It can be seen then that 
for the case of the hollow circular cylinder, as for other confi~ 
gurations, the general thermal stress problem is quite straightfor~ 
ward and evolves itself, for the more complex sonditions of heat 
flow, into that of the determination of temperature distributions 
as a function of time and distance. 

Consider now the more general problem of the thermal stress 


distribution in a cylinder subjected to the flow of a mediua 
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through the cylinder different in temperature to that of the 
cylinder. In this case the inner and outer surface temperatures 
will vary with time and in general will net be linear. The sclu« 
tion of such a problem alee invelves the consideration of a surface 
heat transfer coefficient which itself may have 4 magnitude which 
varies with temperature and consequently with time. Analytic selue 
tions for the transient temperature distribution is a cumbersome 
double series of Bessel functions which does not readily lend itself 
te numerical results. However, analytic solutions are available 
fer unidirectional heat flow in the semi~infinite solid in terns 
of exponentials and the prebability function, which are readily 
adopted to mumerical results and can be aoplied to the hollow 
eireular cylinder within certain limitations. Such solutions are 
valid for infinitely long thin eylinders and give a close approxi- 
mation to the temperature distritution for thick cylinders, in which 
the wall thickness is not large compared to the redius. The use~ 
fulness and practicability of such a method to determine transient 
thermal stresses in hollew cylinders when subjected to a heating 
medium is the subject of this paper. 

In the development of the problem, two separate methods of 
attack were considered and their results compared. First, a numer~ 
ical method for the transient temperature distribution baged on a 


difference equation rather than the applicable differential equation 
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was used to determine the time-temperature gradient. Sesondly, 

@ purely analytic method wes used to determine the time~-tempera~ 
ture relations and in both instances a numerical integrating procese 
was used to evaluate the integrals in the stress equations. It 
might be noted at this point that the atress equations will give 
exact results, within the accuracy of the numerical integration, 

if the time-temperature relations thus determined are the true 
temperature gradients for the heat flow considered. 

It is the intent of this paper to indicate a practical method 
of solution for such thermal stress problems rather than to obtain 
accurate numerical results for any specific problem. However, in 
order to more clearly indicate the feasibility of this approach, 
® simple numerical problem has been carried through, using average 
values of thermal and elastic properties with slide-rule accuracy. 
The solutions of more complex problems are indicated as the method 
is developed. The problem chosen for ite simplicity is that of a 
hypothetical ten~inch I. D. by twelve~ineh 0. D. steam line, ini~ 
tially at 100% suddenly subjected to a saturated steam flow having 
& tenperature of 500%. It is assumed that the steam starts to 
condense immediately and remains at the saturation temperature 
throughout the transient heat transfer process. The surface heat 
transfer coefficient was taken as an average value over the inner 


surface of the pipe, considering film type condensation, and assuming 


. pd we ¢ Pt of 1 ssapegpnae i ef Paes 2 at oa a pore yrrete 


~arequei-ontt od? entmoted of Beau wey bodien ofttiana vi 38 

wetgiereemes eeotuee aoden bee Weed aff nag wht ot 
meee ee Cactseain 4 eaonadae! tod a Son ecoldater ent 

ping ac? . 3 ihe og doe Ay fapeoe, of Rag ai? aely cue 


et saroitaupe toerte add at elaryedat ed etactave of boar ene 
tee itis © te widiwidense off eeetomkt ante anid @ dace te. rival 
ovts Liv ene taupe auotse oid dad? tatog atdd ta Seton od 
tpigv +h « awed ten Sietdh datew drobelitese tote ree i 
soo ktarrg0 3m} ‘aobyemma odd 20 yYowtooss edd mhithe ,edleaet 
ing airyle omit Sile ylanerpeamee hes ili deal ae er nnd 
earit “a. ets Senduseted exit eyes 
eoerodiem a af neltoitebet) enadanemeet ¢ aaa d hanpeal ow oar 


| aberobtaaoe welt aed auld 10> 8 


eoee & (ii ies ton Sag Antty paodioney tegeme ty - A howe wy 


hodton Lao isang # etsotbe? of roqae ett? teodnt edt a2 dt 
ofdelicws oy umidnine @hitiaen ,serewes .stieeet Leslee Of 19 


aiaséo of ant? vodter omaidorg aeetde Samted? dove v0 solduios <3 
i ne 


wave? »s mwiinitiwieldseq wit ai volt dead Jewatteer tila 


al toreved “wo sdoxe ot thenge ‘se tet efiveot Se RO a i 
Gahieet GAs Relite ,apira wt wi hfLidedece Gn? oe. ehsiur ayeete { 
sdonotaes elds Ye wetsidiae? as atastint vinaato own o@ whee 
wei | a (hon og ea? Se inns BO DOK at tye awe lees os Sage Oe i 
contra salt geet Bebveae cand aa maldong Laotvonms oteate a 
te eeobeeios ; 0 oar daviotl efadres aiviviv wielilys S ¥ 
“Wpmerned Sleeennte: 0750 aateeayey ahheiin’ Ree telnet We to sonlay 
uw tvebaiioe aid’ geet tietiod tal cot Balee 
boddon ats so couatee on setae ort Sea 
bays etebalreo Mein? 2o% codtedietedh ewmtaeeee ade at eee 
a 1 Mo dass ye. reemgee. rye wet eset yews eh 
awrel ¢ow «ef AveraR es Tine aad 
-int anks note Re 0 doxt-oviews “W a we" paige fast 


¥, 


aur a mPnh we é Cove Se weiiidealésae bea aueeet 


antvad voi? maoze dodwnuiae @ 03 desvelsom yLnebiwe POO aa Late 


fom et fy on ee a es phage ye Lanaeegett 


ot advade anete odt dadt hommes ob #2 POO to eustereyued 8 


sogaq eid Ip soieioe ale of Ss alehiaal 


onaterrngunt Sanbeneies auld ta ontanos baw hrsmercne gol 


sit : mex at Act w 
pee tt Pe {4 


tao soot ott seanoony zo tans? ‘taod ‘tnoteness was swoon 
{tL wee Gre 


1 


moast sit ove suey epeceva ae wa nese: anv Snetoltieue setanand he 
Madd sear ade Ce Boeew feed ca 
spriyinuon bres ciathoneient on le teary eet ont to vonedn 73 Bis 


eS Past hee Reais" Ae he 


it not effected by temperature, “he magnitude of the surface 
heat transfer coefficient was chosen from averace experimental 
values fer film type condensation as found in the literature. The 
reason for these basic assumptions and their validity are further 


@iecussed in the body of the text. 
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GENERALIZATION OF PROBLEM AND BASIC ASSUMPTIONS 


The general differential equation of transient heat flow is 
obtained from the familiar basic Fourier equation for the conduce 
tion of heat (4 = ~Aeagt at) and a heat balance. Disregarding 
heat sources or sinks and equating the heat gain to the difference 
between the heat entering and heat leaving an elemental volume, 
the general differential equation expressed in rectangular coerdine 


ates ia: 
te 5 T 
3 « &e | 3: (ba3t) +B ley 3 ) §; (235) | 


If k, the thermal conductivity, is considered inverable with 
temperature change and the material is homogeneous and isotropic, 
k will be the same in all directions and constant. 

Then; 
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In general, engineering materials can be considered isotropic 
and homogeneous and this is an ordinary assumption made in engi- 
neerins caleulations. However, the assumption that KE is constant 
may not be valid for all materials within the range of accuracy 


desired in the resulting calculations. The order of magnitude of 
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the variation of k with temperature may be seen by citing mild 
steel ss an example, which hae a thermal conductivity of approxi~ 
mately 26 Bfu/hr ft °F at 32% and 21 BTU/hr ft °F at 1100%. 

The desired tine-temparature<distance relations fer the heat- 
ing or cooling of various shapes are obtained by the integration 
of the general conduction equation, substituting in the necessary 
Doundary conditions for the case in question. One of the simplest 
cases is the unidirectional heet conduction in a plete of semi- 
infinite thickness, with the frees surface at x#@. If the 
term k/ey@ ie defined as the thermal diffusivity and indicated by 
the aymbol K the general conduction equation for unidirectional 


heat flew reduces to: 
oT. k (oT ) 
ot ox? 
fo solve thie equation, it is necessary to determine the 
formulae which will express the initial and boundary conditions 
which the temperature satisfies. If a semi-infinite solid with 
ite plene face at x = 5, initially at the constant temperature of 
tere and suddenly subjected to a heating medium, which has allied 
with it a surface heat transfer coefficient EH, the initial and 
boundary conditions ere: 
Initial condition T=0 WHew t=o 
Boundary condition — at + 4 (15 -Ta) =o. AT x=O 


oqued Mitky x Ye uetsernav ety 


_ A eee seen ote va 
<haamegen to; webeithnntnte: teanedpeentaranaias aiaanadaadaiaas 


1 ES SNS AF Bee, CE Se NP SARE Mee 


aT 


~tasd ont mo exottafor nonstath-orrtet Aetteed edt. 


not gengetat odd we beatatée exe mee i ad 2, wh loge ae yak 


Ay] ae 4 oe ls 


Tine weoen eda at aattod htudire os damp e patteadnns, Lerrenes, ont | 


feeignte sft Yo eno .cotteeny mt gan edt wor enotd tices yamhouod 


~tmoe % etaly a at ae Lhambaep teod | faaet doer thtaw ond et noe 
nat ot +O @ X #8 ooeteye eott sdé déiw ,erendtoidé edkn tink 


Wd Sosaokbat ham yivievttts Lamreds eds "! bent teb at Jqo\x mae 
fanoltyextstne: ‘0% cet fane mottoshers Leveney we 2 Setles on? 
tot ennnter wolt tae 

6)” W dP FG S. alt oak : 
ARE fe a ok ox erapgah 
ot sxtmsotes ot yrsaussen Bt #¥ ,xeteanos oft! e¥foe of 
anoltthaes yrahawed bas fstetat odd meotqxe £ftw dotdy eaLeate% 
tiv bilo od tn} tidaned, s ot ‘Wéetbabbes orwtaneqnos edt dgtedw 
Yo etetareqme? ¢aatenon edt ta yLlatiiat ,4 @ % de enat enata att 
bollia wat doltdy ,arthea antieod % of bodooldve yirebhira baa over 


hos Ishvtnt edd. JX $ootetTiepo setwagx? teed. sondage. 2 92 Aide, 


ne EO TD 
(ond ‘aaah, O=T aolipbnoo Letthad 
Oex TA o.2(eT at) & H+ Tf - antennae, 


er defining H/k as h the boundary condition is: 
= a + RARE TH) eo At x=0 


The solution of the differential equation for unidirectional 
heat conduction with these initial and boundary conditions will give 
the temperature as a function of time and distance for the semi- 
infinite solid. Such a solution will alse be valid for a plate 
ef finite thickness having a large ratio of surface to thickness, 
if it is heated, as before, at the face x = @ and the other face 
is perfectly insulated. Po use such a solution for a hollow cir- 
eular cylinder, externally insulated, and subjected te the flow 
of a heating medium 1+ must be aseumed that the curved section 
of the cylinder can under certain circumstances be considered as 
a flat insulated plate. This assumption will necessarily introduce 
an error in the transient temperature distribution thus determined 
for the cylinder and will be a function of the ratio of thickness 
to redius. Thie fact must be considered in the calculation of 
stress distributions for stresses thus determined are only as reli- 
able as the temperature distributions used in their caleulation. 

To verify the validity of this assumption the applicable 
difference equations may be set up for the plane plate and the 
hollew circular cylinder and the results compared. As explained 
under the heading "Numerical Solution of Temperature Dictribution®, 
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the numerical method converges rapidly to the analytic solution. 
his being the case, 8 comparison of numerical methods should alse 
justify the use of the analytic solution for the plane semi-infinite 
solid for stress calculations for the hollow cireular eylinders 
within the same range of seouracy. 

Sekine a section of the plane solid of unit width, Ax thicks 
and dividing it up ee shown in the figure below, & heat balance 


PIGURE A 


k(T-Th)at + k (T3-Ta) at = Cpe Ax? (Tr - he) 
WHERE: To) = T2 At (t+at) 
Oe (T-h) + (T3-T2) =™ (T'-T2) 
WrHeERe: M = Cela’ 


THEN: GB=eh+th-2zaht+ls 0) 


wa | 
=H + T-2H+ Whee M=4 
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A eimilear heat balence for the circular section where the 


area between sections progressively increases gives: 


K(T,-h)at + & (0 + OXA6)(T3-T2) at = Ce Pax? (\ + 9408)(Toy 7% 
WHERE! Tee = Te Apter (t+at) in CyLinDER 
OR (t-h) + (1.0 +Axa6)(TR-T) =‘ (i+ axae) (Ter -T) 
Z 


SowinG For Ter : LETTING m= 4. 


Try = + fe - Tz (2+4x% 40) + Ts (1+ avao) @® 


mi + AXD0) 
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It can be seen that equations 1 and 2 are identical in form, 
the only difference being the addition of the term 4xoe for the 
eireular cylinder. Por the ordinary problem, in order te satisfy 
the basic assumptions of such a numerical solution, Ax mist be 
chosen as & low order term. This means that the radius of the 


cylinder will be large as compared te Ax and thus Ae must aise 
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be & low order term. The magnitude af these values may readily 
be seen by taking as an example a hellew cylinder of five inch 
inside radius and a ox = .01 ft (Hote that this is the largest 
value ofax chosen for numerical calculation) then: 


A® = TanaAao = x2 = .024% RAD. 


and Ax4e = .00024 (Whieh is essentially a 
second order term). 
Then for the use of a modulus M 2 4; 


i] 
‘Tee = Te. 4 [t - (2.00024)T + (1.00024) Ta (©) 
4.00048 


A comparison of equations 3 and 1 shews that thig point, 
nemely at a radius of five inches, Ty and t. are practically 
identical. However, in the mumerical solution for temperature 
distribution, the temperature at any point is used to determine 
the temperature at the next adjacent point, ete., through the 
solid. This means that the error in using equation 1 instead of 
equation 2 will be an accumulative error and places a restriction 
upon the thickness of the cylinder for a given error. This fact 
perhaps can be seen more clearly by noting from Figure A that; 

AXn =A¥+ HaXo6e 4 
Where: n is the number of steps te be taken 


for the given thickness shosen. 
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It should alee be noted that the error is not a function 
of Ax alone but a function of the product of ax and AO, 
which essentially means that for a given radius cylinder and wall 
thiekness the error involved will be the same regardless of the 
choice of ax. With these facta in mind, it ie apparent that the 
limiting factor is the ratio of wall thickness to radius regardless 
ef their absolute magnitudes. 

Another pertinent observation ise that at the first instant 
of heating the two equations almest coincide and the errer in 
temperature distributions increases with each step. As the tempera~ 
ture distributions — Sens Sy enter the stress equations 
in the terns ft rdy and Trolr the true criterion for 
the error nneee in etress ciliates would be to plot tempera- 
ture times radius vs. radius curves for both methods and sompare 
their areas. 

It is the author's opinion that the error involved due to this 
basic asgunotion will be of the same order of magnitude or leas 
then the errors inherently involved in thermal and elastic proper= 
ties, where ratios of wall thickness to radiue found in ordinary 


ateam lines are to be used. 
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NUMERICAL SOLUTION FOR TEMPERATUNE DISTRIBUTION 


Unidirectional heat flow will occur in a flat slab at 4 
distance from the edges and also in an insulated red of wiiforn 
section. As far as temperature distribution is concerned, the size 
end shape of the rod is immaterial and is only of consequence in 
Galeuleting total heat flew. Consider such o rod of unit width, 

a thickness ax and divided up inte sections Ax long as indicated 


in the figure below! 


oP 


If now a finite surface heat transfer coefficient H is con= 
gidered tetween the surface (0) and the ambient fluid temperature 


(a) then for the point (0) a heat balance gives: 


Hax(Ta-Te)at = - ke (T-Tejat + C@ax® (Te-Te) 


Where 1p! =o Ar (t+at) 
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The finite difference equation thus obtained is based 
upon the following assumptions: 

a That At be taken emali enough so that there is a nege 
ligible error in heat flow calculated by the tern 
(4p )/ox. 

b. That At de taken small enough so that there is a neg» 
ligible error in neglecting any tut adjacent areas 
(ax by ox). 

ec. Thatax be taken amall enough so that the tempersture taken 


at the ceniter of the element can be used to compute the 


heat capacity. 
Define: Mc GoP ax 2. ex? 
RAL K at 
N= Har 
i 


Using these relations: 

N (Ta-Te) +(Ti-Te) = 4 (Te-Te) 
Solving for TA 

tis 2N To +(\- ZNAZ) Te + 21 


This equation may be written in terms of "influence factors" 


ae follows: ‘ e 
Mee Paola + Foo To + Fo Th 
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The subscripts indicate in order, the influence that the 
present temperature has on the point whose temperature at (t+ At) 
is to be found for the period of time At. In other words, if the 
new temperature st point (0) is to be found after the time At, 
Fas indi¢ates the influence of point (a) on the temperature at 
point (0) during the time interval At. 

if fj is to depend upon %, in a positive sense, then: 


Fue +1~ 2Na2 26 
nal 


OR M-2N+2 26 
MZ 2N+2 

Since ail of the quantities entering into HW are inherently 
positive, it is apparent thet M must be larger than two , for the 
case where a surface transfer coefficient is involved. Depending 
upon the magnitude of the surface transfer coefficient and the ther- 
mal constants of the material, this may place # limitation upon a 
convenient averaging process. The value of M=4 was found suitable 
for numerical temperature calculations for steel. 

The general method is to use the weighting process, as described 
shove, for surface temperature calculations and the simpler averae 
ging process for interior temperature determinations. It is con~ 


venient to choose a scale euch that fractions can tbe neglected. 
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If the fraction 1/2 was consistently added or dropped, an error 
would tend te aceumilate. ‘Therefore, a rule of dropping the 1/2 
efter even digits and adding 1/2 efter odd digits can be adopted, 
ast 1/20, 3/22, 5/22, 7/24, ste. 

The choice of the modulue M is discussed in reference (5). 
The larger the modulus M the sooner the numerical results will 
converge to the analytic solution (i.e. in less steps) end the 
greater the number of calculations necessary to cover the same time 
period. After due consideration of the work involved, @ modulus 
M=k was chosen to obtain rapid convergence to the analytic solution 
a@uring the early time intervals. 

Now consider the heating of a one-inch thick plate insulated 
on one surface and subjected to a heating medium on the other sure 
face. Yor simplified calculation, consider E and H as constant. 
By choosing several values of ox (and the compatible at) for a 
modulus M=4, surface temperature vs. time curves can be plotted to 
determine the validity of the choice of the magnitude of az. As 
the accuracy of the solution is based on the assumption that ax 
be emall, such a plot should indicate any variation in the results, 
due te an increase in the magnitude of ax. 

For ordinary Steal: 
K=.50 £7 /ne 
k = 25 Bru/hr, ft, 
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Assumes H= 1000 Bru/hr, ft*, °F (A fair approxina- 
tion for film~type condensing steam -~ see Appendix 1) 

Sample determination of My Fags Foos Fo, and at. 

iat AX = .01 FT. ay = 


=12 x.0o1 = .i2 INCHES 
Then: N = 


Hax = \000xX.01 — .40 
Ae 2.5 
M= ax On At= ae WHere M=4 
Kat 
at = (0% = 0.5 x10 Hes. 
4x.5 
At = 0.5x10°*x 3600 = .18 SEC. 
Fie — 2N -z 2% 0.4 2 0.2 
~ 4 
Fea is N42) = (\-9.8+2-0) =0.3 
~ 4 
Fie = Zz = 2 = os 
™m 4 


Similarly for the other values of Ax: 
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The equations for the surface temperature become; 


a? B= of + 36% 


- 0.5 T, 
(2) Tt = Git © 6-35 Tos + 0.57, 
= alcta + een «+ cer 
) te = 0.05TA + 045T. + 0.57, 


On set of calculations will be given at this point te indi-~ 


eate the form and procedure. The results only for other numerical 


ealoulations are included in Appendix IT and plotted on the enclosed 


curve sheets. 
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ANALYTIC SOLUTION FOR TEMPERATURE DISTRIBUTION 


fo determine the time-temperature-distance relations in the 
eemi~infinite solid for the case where the initial temperature is 
sero and the surface temperature is a function of time $(t), the 


following equation is valid: 


i, 
2 
=A 
yap b(t -—x je daw ® 
af ( 4K 2 
X/2VKt 
Also for the case of the semi-infinite solid, initially at 


sero temperature, which is heated at the surface x=0 by radiation 


from a medium at constant temperature Ty, the applicable equation 


ie: 
hu th? Kt 
ie sents Ms ote eric 4 M4 be OY (Cz) 
"Ta 2VKt 2V«t 


Where erfc is defined ae one minus the probability function 


or error integral. 


evfe x =(1- ae fe sn) 
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It can readily be shown by direct substitution into the 
differential equation of unidirectional heat flow that these 
eqaations do satisfy the differential equation and initial and 
boundary conditions. 

Zquation (1) was civen to emphasize the fact that whenever 
the surface temperature of the dildbinsittnans $0 Meth as a function 
of time $(t), and this function is continuous and can be integrated 
ever the time interval (in equation (1)) to be considered, the 
resulting réneneipanntuite<thuteene relations for the saclid are 
then determined. hie opens the avenue for the solution of a number 
of interesting problems. pe noted before, a numerical solution 
for surface temperature is readily adopted to thermal properties 
of the material which vary with temperature, as well as a surface 
heat transfer that te a function of temperature. This method can 
also be extended to the case where the temperature of the heating 
medium is some known function of time. However, it should be noted 
that such a solution for surface temperature, under the above stated 
conditions may become very tedious, as the influence coefficients 
may have to be calculated for each step of the process as well as 
choosing a At that ie small. Although the numerical solution 
iteelf, if carried through the solid to the desired depth, will 
give the time-temperature-distance relations directly, the desired 
surface temperature $({t) may be determined over a shorter interval 


with the resulting saving of much numerical work. 
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Figures (c) and (D) show the determination of the surface 
temporature relations o(t) jor a short interval of tims. Figure 
(D) shows that the surface temperature ve. time curves very closely 
epproximate a straight line on & log-log plot, and can readily be 
pat inte an analytic expression. 

Equation (2) has been used to determine the time~tomperature- 
distance relations in a perfectly-inaulated, one-inch-thick plane 
#olid of large surface area, The results of this can be used as 
a close approximation of the true transient temperature gradient 
in an insulated 1l0-ineh inside diameter hollow cylinder of one~inch 
wall thickness. 

Yor the ultimate caleuletion of the stress distribution, it 
ie convenient to plot temerature-distance relations for variqus 
intervals of time, in order te perform a graphic or numerical inte- 
gration of the integrals in the strese equation. Figure (8) shows 
& plot of these relations and shows how Glosely the numerical method 
parallels the analytic solution. Figure (F) gives the equivalent 
temperature times distance vs. distance relations, and indicates a 
arephic integrating process. It might be preferable and possibly 
@ little more accurate te use a numerical integrating scheme sach 
as the trapezoidal or Simmson's rule. 

Wunerical regulte are not as tedious as a first inspection of 


equation (2) might indicate. Tabalated values of the error integral 
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are available and many cf the parameters are constant for one 
instant of time. A simple tabulation of the quantities given 
below will materially facilitate calculation, all other para- 
meters being constant for a given instant of time. 


(1) X% iwFr 


(2) x 
VKt 
(3) ev fc er 
(4) evfc (ee + h VKt} 
(5) nx 


(6) hx + h?kt 
(7) nxt nekt 


(8) sf 


The results of calculations made for the semi~infinite solid 
for a penetration of one inch are given in appendix to correlate 


the data used to plot curves. 
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DETREMINATION OF STHESS DISTRIBUTION 


If the transient temperature distributions determined by either 
numerical or analytic methods for the semi-infinite solid under the 
conditions specified are considered good approximations of true 
temperature distritmtionse in the hollow cylinder, caleulation of 
stress distributions are etraight forward. At any instant of time 
the stress becomes a function of the thermal and elastic properties 
of the annentes for the ceometric configuration considered, and can 
be readily caloulated. The following thermal stress equetions for 


the finite hollow cylinder are given 4 Timoshenko: 


= Se. JP eee ps ss pd 
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These equations are based on thermal and elastic properties 


i 
1% 
xT 


that do not vary with temperature. Within this limitation they 
are exact for an infinitely long cylinder of any wall thickness 
and any given temperature distribution eo long as it ia in the radial 


4irection only. The principal planes of stress are in the radial, 
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exial, and tangential directions with normal stress acting on 
each and stress acting along a plane perpendicular to this axis 
are independent of the pesition of this plane. 

In order that the solution be exact when a cylinder of finite 
length be chesen, boundary conditions require that external forces 
equal in value and distribution to the stresses in the axial direction 
act on the ends. In the usual case such forces will not be present, 
so that the stresses near the enda will be very different than those 
calculated by this method. ‘This effect diminishes as the distance 
from the ends is increased and practically disappears at a point 
sufficiently far away. Writers differ as to the distance from the 
ends where this point is located and it actually depends upon the 
ratio of the thickness to cylinder diameter, For the purpose of 
our problem, this influence will be considered negligible. 

I% ie interesting to note that in the case where the eurface 
of & solid, initially at uniform temperature, is euddenly raised to 
a given temperature and thus maintained, the maximum etress occurs 
at the surface at that instant. In fact, the maxinum stress (i.e. 
tangential and axiel stresses which are equal at the surface) will 
always occur at the heated surface for a solid initially at uniform 
temperature. The stresses tims produced are equivalent to the stress 
produced in a thin plate for the same temperature rise in which & 


chance of area is prevented. This stress is only momentary and 
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rapidly drops as the increase of temperature spreads inwerd and 
allows an expansion of the outer layer. Most published solutions 
have this common boundary condition, namely, an instantaneous 
surface heat source. 

For the case of radiation inte a solid from a fluid, where 
the surface temperature is a function of time and equal to the 
initial solid temperature at the first instant of heating, the 
conditions are vastly different. The maximum surface stresses 
will not occur at the first instant of heating tut at some later 
time, and depending upon the conditions of heat flow may even occur 
at the steady state condition. Por the case of the sudden rise of 
the surface temperature of s solid, initially at constant tempera- 
ture and ultimately coming to a steady state condition, it is essen~ 
tially going from an unstressed condition to a stressed condition. 
In the problem chosen, however, the cylinder at an initial uniform 
temperature with the outer surface perfectly insulated will, upon 
heating of the inner surface, be taken from one unstressed condi- 
tion te another as the cylinder finally tinshen the tempsrature 
of the heating medium. It is reasonable to essume then that maxi~ 
mum atress conditions will occur at some intermediate time. 

Calenlated values of inner surface strees are tabulated in 


Appendix II for a one-inch thick, ten-ineh I. D. circular steel 
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cylinder along with the stress distribution ecross ite section at 
the time of maximum stress, The cylinder being insulated on its 
outer surface, initially at 100°F and suddenly subjected to a satur- 
ated steam flow which is maintained at 500°F. The results of these 
ealeulations are graphically represented by figures (£), (-) and (G) 
and will be further diseuseed in the results and conclusions. 
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APPENDIX i 


Accurate thermal stress determinations in eteam lines are 
usually linked with a knowledge of the surface heat transfer co~ 
efficient. The matter of the physical phenomena involved as well 
as the magnitude of such/ coefficieats is the subject of a good 
deal of conjecture. This coefficient is defined by the heat flew 
equations 

ge HA(T,g-Ta) 
Wheres q ~ Heat Flow 
H - Surface transfer coefficient 
A + Transfer area 
Tg- Fluid temperature 
T—g- Solid surface temperature 

Values of H determined by experimental means are in general 
more reliable than these determined analytically, and particularly 
fer condensing stean. 

Phere are two types of condensation, namely, dropwise end 
film condensation. Zither or a combination of the two types may 
exist, depending largely woon the cleanlinesa of the fluid, cleane 
liness of the condensing surface and the surface conditions of the 
material concerned in the condensation process. In any case, Leat 
transfer coefficients for the dropwise process are in general many 
times higher than for film condensation, but are more difficult 


te determine, 


GW ISS VIR GantCs ae Mer PAE ade ait apbe prin fa Risto 
ZS . 
ad! ne Satatweed? gadat ete #2 . 9 eeeaa aetna te wate er 
wise & af heiealdor elohine few Sos ip liad i far went neg 
oa pain maeta at anoltantnxe2es snouts Laerredé 
rade SS atte wet 0 it de betdetelbew of unkde walt ands bebe 
“90 ae tener tear sae ture ad? to esbelverd * détw bedats Yi lawew 
ey eek ee Le ioome ert ©) party, sent wrens 


tte #6 sevdowat acsaonase tastes oat %» atta oe ae ed LLie im eee 
boon @ % dootdwe ead? at siaetottieos ‘siowe Yo shatiogam ad? se ee 
wolt seed ed? yd beatteb of taeiolYtees eldl .eiésetace to iaeb a ie 
ER oe oe 
welt sack Pp tered il : 2 a 
tmetottiee0 retensts osu + H to 
asta velanat = A | : 
sustsregm? som tue SELok =gf | 
latenes si ete eseon Latnontroqae Yd Soateresed K %0 seutat 
Ultelootiaeg bas sClilealiylens Seaierreteb eased? amlé eldaifer exom ‘ 
be selvqoub ‘andi juadhcabaiabeey nah arate hineal 
We sogtd ow! edt Yo sottesidaes 2 to tedt2® .meideenebaoo mfit 
~aiels ,Sief? sit Yo saentinasio edt sogr ylogtel gethaewet ,tebxe 
add to enotithnen soatume edd bas esa trme saleneheoe end Yo woot! 
jaed ,seao Yas cl ,aseceTg molieaneheos oft wi Besrvennen tabvetan 
Visa latesag af @45 aaesotg eulwgorh oft so? atmelolYisee +o tener? 


HlosiYths oxom ois dud ,soldaegobace wilt set madd xesintd eewhs 
-suilateteh of 


Authorea differ widely as to the magnitude of the difference between 
the two types, giving coefficients from 4 te 20 or more times larger 
for the more efficient dropwise process. It bas been deterained 
that in coneral clean steam with or without non-condensibdle gases 
always condenses in a film on clean surfaces, whether rough or 
smooth. 

Experimentally determined values of heat transfer coefficients 
are given in terme of an averace value over &® prescribed area, 
discounting the fact that the value ef the coefficient may vary 
considerably ever the area depending upon the film thickness. Wusselt 
in an analytic investigation of heat transfer coeff ietente on eyl- 
indrieal surfaces found the film thickness and consequeatly the 
transfer coefficient to vary around & horizontal tube, His analysis 
shows the minimum film thickness and consequent maximum transfer 
eeefficient te exist a few degrees on either side of the top of a 
horisontel tube. His investigation is equelly valid fer condensa- 
tion inside or ontside of a tube. By integration around the tube 
an average coefficient was determined. Experimental determinations 
have shown the average ceefficient thus determined toa be somewhat 
in error and correction factor curves are available for use with 
Musselt’s equation to give representative values of the surface hoat 


trensfer ceefficient for condeneing atean on tubes. 
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If now dropwise condensation were to be considered, the 
stress problem immediately becomes very involved for the lecal 
transfer coefficient would vary appreciably frem droplet te un- 
wetted surface end the heat flow would ne longer be unidirectional. 
Por auch a problem an entirely different approach would have to be 
considered. The value of H-1000 Bru/hr ft® OF used in the mmerical 
ealeulations is an arbitrary value chosen from the references for 


film-type condensation of stenm on a horizontal tube. 
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RESULTS AND CONCTIISIONS 

It is apparent from the magnitude of the inner surface strese 
determined, for the wumerical problem considered, that thermal 
atresses wader Severe Gonditions of repeated heating cyeles could 
eontribute to ultimate rupture of @ line under pressure. The tan~ 
sential atresa on the inner surface due te preagure in the Lins 
will be tenelle and to some degree domtemes thermal stresses. 
However, for the l0-inch I, D. by 12-inch 0. BD. line considered, 
the stress due to a pressure of 515 #/in* in the line would be Bim 
proximately dy = $00 #/in® in compression and ve = 1,640 #/in® ia 
tension, These atresees are extremely omall in magnitude as come 
pared to the thermal stresses determined. 

It should be remembered that this paper was motivated by a desire 
te formilate a practical approach to the actual theraal stress 
problems in hollow cireuler eylinders which are extremely diffieult 
if not impossible by purely analytic appreach. Necessarily, many 
aseumptione were made, and to some degree they have been justified 
in the text of the paper. It 1s obvious then that the reaults of 
the numerical calculations should be substantiated by an exper iareric 
tal investigation before they are accepted as accurate enough for 
practical use. 

4s little data is available on the magnitude er variation of 


surface transfer coefficients, fer condensing stean in hollow 
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ciroular cylinders under high temperatures and pressures, the 
value chosen could be considerably in error, a8 it is a value 
extrapolated from experimental results for condenser tubes under 
lower temperatures and pressures. It is conceivable that a 
gensideretion of the change of the magnitude of the surface trans 
fer coefficient with temperature difference across the film as 
well as the change in thermal properties with solid temperature 
would waterially affect the thermal stresses determined. The 
simplifying assumption that the line was perfectly ineulated would 
alee constitute some error in stress calewlations beyond the tine 
at which the outer surface temperature begins te rise. Although 
under steady state conditions ineulations are availeble that are 
extremely effective, under transient conditions, temperature ne 
tribations will necessarily be effected by the heat capacity of 
the inwuletion, This condition was net considered in the numerie 
eal problem which injects some error in the resulting stress deter- 
minations for times beyond approximately five seconds. Perhaps, 
with ome of the qualifying assumptions in mind, this paper will 
stimulate further investigations along the same lines teward a 
usable practical gelution for thermal stresses in steam lines. 


os 68 


Quer e seth Sa ‘s oA ae NP cy 
8 tat exaurivoaes 0 Mh. 


Weer 5 SS Rey GT ME phe es CRN ‘ 
~saats soatuve odd to eb 
by xem Bees * Kd Bai ca eM pe eee 


ponte SPAS, eit Me, Re Tie 
hiuow boratuent “toon a eats mt ae pireaaiiongtatec 
only ot boxed enol tatentae svente at fh soeee amps, ated ttonee, ope 
Rare Merde 2: snetaem metre, ot, Sadey Ae 


‘ire 


Xe 


gf ) 5, F & eh Aes so ceeraters’ 


“sath erode: > mnie, dentenert, abey santero Sempre 
be ban So 2 oe 8 po A ¥ Seteete | "., ee SR fhe ano lindbsd 
~bsenar ae - bexabtexoo we ms ‘ena ome sda Laren wate 


ste nnn anton dt ok ore mo outed dnt ate, Sm 
segaiise% vabacees ovd? vletenduonees haoyed eenit wot enotéentm 


PMOL He PS bee rey ERTS wegen f 


{et sega ais sce af anton antec ot % ome iw 


me aa. = ev ataey, wre eh oi aes @ 


a bumeod secs oem sd yeotn anoteapttacvat mut 


sf, 


Ly Feed gaat td 
sooclt mages at asnowtée fama wei aoidalen Laokinere ofan 

Paton Laces Gate 

Pa MS | me a 4 er Boxe orld Cyt 


: 
, 
| 


RECOMMENDATIONS FOR FURTHER STUDY 


Further study of the practical problem of thermal atresses 
in hollow circular eylinders during actual transient boundary cone 
ditions falls inte the several brosd categories listed below: 

(a) The experimental determination of the true boundary 

conditions thet may exist in the actual problem during the 

transient state, such ae those conditions that may exist 
when a steam line is suddenly subjected to a flew of high 
temperature steam. 

(>) An experimental verification of the numerical example 

given in the text as to the time at which maxinmm stress 

ooeurs, ite magnitude, as well as the stress distribution 
across the section at time of maximum stress. 

(e) An experimental study of the Magnitude of the surface 

transfer cocfficient under high temperatures and pressures 

in holiew eclreular cylinders. How the value of such a cox 
efficient varies with proesure, temperature difference 
acrose the condensate film, and aneular position in a horie 
sontal cylinder. 

(4) A further justification of the use of temperature dis- 

tributions as determined for unidirectional heat Plow in 

the stress caleuletions for hollew circular cylinders. 
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There are two imnedinte avenuca of approach to the last 
dategory, first, by direct comparison of the temperature distri- 
butions as detersined by purely mumeriecal methods for the plane 
plate and hollow circuler cylinder, and, secondly, such a tempera- 
ture comparison could be made by the use of a conformal mapping 
process from a plane finite section to a cireuler section. A plot 
could then be made of the error in stress calculations as a funce 
tion of the ratio of wall thickness to radius, which would verify 
the usefulness of the approximate method for e problem of given 


geometry. 
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APPENDIX IT 


Tabulated values used to plet curves 
Figure (A) « (Numerical method ~ Surface temp. ) 
Te! = Ta + (4 - NE) To +4 (Nore: m=4) 


To =100, Ta = S00, H= 1000, K =.50 


Ax = .12 inches x = .09 inches Ox = .06 inches ax = 03 inches 
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226 lg 192 .270 150 0.06750 
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220 060, 808 203.360 
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Figure (A) = (Analytic method - Surface temp.) 
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To = 100%, tT. = 500%, h = 40, K = .50 
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143.0 0.05 
159.6 0.1 
177.6 0.2 
192.4 0.3 
204.3 0.4 
215.2 0.5 
221.7 6.6 
229.7 0.7 
236.4% 9.8 
242.0 0.9 
247.0 1.0 


Figure (6) ~ (Numerical method ~« Solid temp.) 


Surface temp. -- To = O2Ta + 6.3T. + O5T 
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Figure (@©) « Analytic Method ~ Solid Temp. 


f, = 100%, f= 500°F, H = 1000, k = 25, h > 40, K = 50 


t+ = 0.1 See | $§ 3 0.2 See t+ = 0.5 See $ = 1.0 See 
Radius Radius Radius rete 
T-% in * t- in* t~% in pe “Fy in® 
158 0 45.0 174.8 5.00 214.4 5.00 247.5 $000 
138.9 45.02 146.2 5.03 188.8 5.03 208.7 5.05 
124.6 5.04 130.7 5206 166.2 5.06 176.8 5.10 
1M.1 865.06 116.4 5.09 148.0 4,09 152.0 5.15 
107.6 5.08 108.7 5.12 132.5 5.12 134.0 5.20 
103.8 5,10 104.0 5.15 122.9 5.45 120.5 5.25 
101.8 45.12 101.7 5.18 115.2 5.18 112.1 5.30 
100.8 5.14 100.7 5.22 109.6 4.22 106.8 5.35 
100.26 5.16 100.3 5.24 105.8 5.2% 103.6 5.40 
103.5 5.27 101.7 5045 
101.9 5.30 101.3 5.50 
101.0 545.33 
100.5 5.36 
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Figares (C) and (D) ~ (Mumerieal Method ~ Surface Temp.) 


To’ = NT + (4 — Mt! \Te + oh (Note mM=4) 


Ax =.09,) at= Ot, Te =100, Ta=500, K=.50 


Time in 

See. 3.2 900 #51000 B= 1100 # = 1200 
e 100 100 100 100 
2101 154 160 166 172 
2202 174 181 188 195 
° Ne 188 196 204 211 
e 199 208 216 224 
2505 208 217 226 235 
° 06 216 226 235 2hyly 
2707 223 233 252 
2808 yt UR 249 259 
«909 235 246 255 265 
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Ficure (E) ~ (Anmlytic solution « Temperature times radius 


ve. radius) 

_ t = 1.0 See t = 3.0 See &t = 5.0 Seo + = 7.0 See t+ = 10.0 See 
Inches Oe Oe BRR: linen Raion 
$00 738.0 10338 1175 1276 1352 

° 492.9 758 9 1055 1179 
5,20 177.0 520 eh 865 1008 
5530 64.0 337 527 676 821 
$440 19.5 203 $24 675 
5250 ‘ 1 257 385 
$60 61.0 170 287 is 
5.70 30.2 126 198 

. 9.8 58.6 144 260 
5.90 > 31.9 89.0 188 
6.00 1.2 19.8 46. 7 


‘ 


.. 
4 


Ose DOE = 2 098 MT Ft ORF ORR E gee 0 
er 2 ol 4 as eo ; % iy 5 , he \ " 
oy. oe he 8 a"; ra ! iu 


DET ak 


Figure (G) ~ (Analytic Solution « Surface stress vs. tine) 
Surface 
steess $21.0 Seq $23,0See $25.0 Seg 857.0 See 1510.0 See 
Vos/in® -38,950  =49,700 $0,250 =%9,800 47,800 


Figure (£) ~ (Analytic Solution - Stress distribution at time (5.0 sec) 
ef maxisum stress) 
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90 19,210 0.0 19,2160 
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